Nebular Abundance Errors 



Jason Alexander and Bruce Balick 
Astronomy Department, University of Washington, Seattle, WA 98195 
Email: jalex@astro. Washington. edu,balick@astro. washington.edu 

ABSTRACT 

The errors inherent to the use of the standard "ionization correction factor" 
("■^cf") method of calculating nebular conditions and relative abundances of H, 
He, N, O, Ne, S, and Ar in emission line nebulae have been investigated under 
conditions typical for planetary nebulae. The photoionization code Cloudy 
was used to construct a series of model nebulae with properties spanning the 
range typical of PNe. Its radial "profiles" of bright, frequently observed optical 
emission lines were then summed over a variety of "apertures" to generate sets 
of emission line measurements. These resulting line ratios were processed using 
the icF method to "derive" nebular conditions and abundances. We find that 
for lines which are summed over the entire nebula the i^F^derived abundances 
differ from the input abundances by < 5% for He and O up to < 25% or more 
for Ne, S, and Ar. For resolved observations, however, the discrepancies are 
often much larger and are systematically variable with radius. This effect is 
especially pronounced in low-ionization zones where nitrogen and oxygen are 
neutral or once-ionized such as in FLIERs, ansae and ionization fronts. We 
argue that the reports of stellar-enriched N in the FLIERs of several PNe are 
probably specious. 

1. Introduction 

Chemical abundances of gaseous nebulae are among the most frequently used 
astrophysical tools for measuring the rates of heavy element formation, inferring the history 
of star formation, and gauging the effects of internal stellar nucleosynthesis and mixing. In 
this paper we subject the most common method of measuring heavy element abundances, 
the "ionization correction factor" (""icF") method, to a simple but rigorous examination. 

The icF method was initially developed by Torres-Peimbert & Peimbert (1979) as a 
method for converting measured nebular emission line fluxes into chemical abundances. 
The application of the method is straightforward: Using standard diagnostic techniques. 
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the electron density, n, temperature, T, and the ionic abundances (relative to H+) of species 
with visible lines (e.g. He+, He++, N+, 0+, 0++, Ne++, S+, S++, Ar++) are derived from 
observed line ratios. Total chemical abundances are then derived by correcting the ionic 
abundances for unseen ionization stages using icpS- 

The many untestable assumptions used to formulate these correction factors render 
them the weakest link in this methodology; however, the resulting errors have never been 
stringently assessed. Various statistical studies of nebulae have been used to refine the 
various correction factors of the icp method, the latest and most exhaustive study being 
that of Kingsburgh & Barlow (1994) (hereafter, KB). However, these studies all rely on 
observations of nebulae whose assumed chemical properties arc based on icF techniques - 
a somewhat circular process. More recently, Henry et al. (1996), Kwitter & Henry (1996), 
and Kwitter & Henry (1997) have developed a method in which icF abundances are only 
used as a starting point for Cloudy models (see Ferland et al. 1996 or Ferland et al. 1994 
for a description of Cloudy). Making certain assumptions about the nebular structure 
and stellar properties, Henry et al. (1996) vary the input abundances to Cloudy until its 
predicted line ratios match their observed ones. They find that averaged over the nebula, 
icF abundances differ by < 20% from those used in their most successful Cloudy models. 

If this were the entire story then this paper would not have been submitted. However, 
Balick et al. (1994) observed several planetary nebulae (PNe) with I'.'S resolution along a 
slit. Using the icF method they derived radial increases in nitrogen abundances as high as 
a factor of 7, whereas other elemental abundances remain constant. In addition, Guerrero 
et al. (1995) have drawn similar conclusions about nitrogen variations in the bipolar nebula 
Ml-75 while Corradi et al. (1997) found radial increases for hehum, neon, sulfur, and argon 
in their long-sht studies of IC 4406. Based on such abundance variations with radius 
some authors have proposed the recent ejection of high-velocity, element-enriched material 
by the central star. If real, these interpretations suggest some bizarre stellar enrichment 
process for which no theoretical understanding exists. 

Alternately, the icp-derived abundance gradients could be erroneous. Investigating 
this possibility is the goal of the present paper. As we shall see. Cloudy simulations 
show that there are important regions, such as all low-ionization zones, in which the icF 
abundances of some elements can be highly discrepant (a factor of three or more). Other 
total abundances befall similar though less extreme fates. 
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2. Methodology 

We denote elements by standard abbreviations, H, He, C, etc. Unless implied otherwise 
in context, we shall use the following abbreviations for specific emission lines: [Nil] = 
[Nii]A6583A, [Oi] = [Oi]A6300A, [On] = [Oii]A3727A= [Oii]AA(3726 + 3729)A, [Oiii] = 
[Olll]A5007A, [Nelll] = [Nelll]A3869A, [Ann] = [Arlll]A7136A, [Ariv] = [Arlv]A4740A, 
[Arv] = [Arv]A7005A, [Sll] = [Sll]AA(6717 + 6731)A, and [Slll] = [Slll]A6312A. Generally 
speaking, total abundances are relative to H and ionic abundances are relative to H+, so 
that "Ne" means the Ne/H abundance and "0++" is equivalent to 0++/H+, etc. "Profile" 
implies the radial distribution of some quantity, e.g. an emission line profile. The electron 
density, n(S"'"), is determined using the 6717A and 673lA lines of S"*" and the electron 
temperature is determined two ways; T(N"'") using the [Nil] 6583A, 6548A and 5755A lines 
and T{0~^^) using the [Olll] 4959A, 5007A and 4363A lines using standard methods. 

In the present study we used Cloudy (version 90.02) in order to conduct a controlled 
study of icF^derived abundances and their sensitivity to changes in nebular conditions 
such as nebular geometry, stellar temperature, luminosity and emergent spectrum, and 
abundances. This procedure bypasses all observational problems (e.g. reddening corrections 
and fiux calibrations) and does not rely on a statistical or empirical analysis of the 
observations of a heterogeneous sample of observed objects. 

Based on a set of assumables" , namely, the chemical abundances A{X) relative 
to H, along with a nebular density structure, n(H,r), and a standard star model, 
Q^{T^,L^) or Q*{T^,g^), Cloudy was used to compute the profiles of the ionization 
structure, electron density, electron temperature, and the emissivities of bright optical lines. 
Cloudy's emissivity profiles were summed over different "apertures" to synthesize a set of 
^^observables" , i.e. emission lines fluxes in relatively bright nebular lines. 

The original proflles and these observables were normalized to H/3 and then used to 
determine a set of ^^derivables" (denoted with a prime) from ratios of observable lines - 
n'(S"'"), T'(N+), T'(0++), and ionic abundances A'{Xion) - using standard techniques. The 
expressions for this process were obtained by flts to output from the Ionic and Temden 
tasks in the IRAF/STSDASQ program Nebular (Shaw & Dufour 1995). These flts, which 
are generally accurate to 2% or better, are listed in Appendix A. In all cases we used T(N"'") 
and T{0~^~^) for computing the abundances of ions emitted from low- and high-ionization 
regions, respectively. 



^Image Reduction and Analysis Facility is distributed by National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 
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The derived ionic abundances, A'{Xion), were used as input to the icF niethod used by 
KB in order to derive the resultant total abundances A'{X). To this end, the A'{Xion) are 
multiplied by certain icF corrections, icri^), which are based on geometric assumptions 
about the relative volumes occupied by the unseen ionization stages. Expressions for 
the icF correction factors for various atoms are also listed in the Appendix. Finally, the 
"discrepancy ratios", A{X) = A'{X)/A{X) are computed. (Note that our A{X) = l/i{X) 
of Henry et al. 1996). 

As an example of this process, in many nebulae the dominant ionization stages of 
nitrogen arc and N^"''. Although N+ emits several bright optical lines (rendering N'^/H"'" 
readily measurable), does not. Hence we must estimate N by correcting N+ with an 
icF where icF(N) =N/N+, recalling of course that N is not directly measurable. In this 
case, the icp method assumes that both nitrogen and oxygen have the same ionization 
structure; i.e. the fractional ionizations x(N+) = x(0+) and x(N++) = x(0+'*") at all 
radii. Thus ic'F(N) = 0/0'^, where (0/0+) can be derived from observables, as shown 
in Appendix A. Note that no explicit knowledge of the nebular physical conditions or the 
excitation source is required beyond that needed to determine ionic abundances. 

Using this technique, an assessment of the global icF nebular abundances was made by 
integrating the line emissivity profiles computed by Cloudy over four regions or apertures: 
the entire model nebula (truncated where a;(H+) < 30%); a "high-ionization" region (in 
which a;(0++0°) < 50%), a "low-ionization" region (in which a;(0+) > 50%), and a 
"neutral" region (in which x{0°) > 50% and a;(H+) > 30%). The high-ionization region 
is typical of density-bounded nebulae observed with a large aperture. The low-ionization 
and neutral regions are characteristic of low-ionization features, such as ansae and FLIERs, 
commonly found in planetary nebulae. The computations stop at ^^(H"'") = 30% since the 
assumptions of ionization and thermal equilibrium become poor at about that point and 
the use of Cloudy is not appropriate. (Even at x(H+) 50% the detailed results are 
probably uncertain.) 

The procedure was apphed to a "baseline model" in which the total density 
is cm~^ within a radius Ro — 10^^ cm, and 3000 cm~^ beyond it. Abundances 
logA(H:He:C:N:0:Ne:S:Ar) = 0.0 : -1.0 : -3.52 : -4.0 : -3.22 : -3.82 : -4.82 : -5.0 
were used. A star with logL* = 37.4 (erg s~^)and a blackbody spectrum with temperature 
logT^ = 5 (°K) was adopted. This baseline model is representative of typical bright galactic 
planetary nebulae such as NGC 3242 and 7662. 

Several input parameters were then varied in order to ascertain how the derived 
abundances A'{X) and icF discrepancies A{X) are affected. The parameter variation 
studies are intended to explore reasonable ranges of structural and chemical conditions. 
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These parameter variation sets, A, are summarized in the left column of Table |l|. 

The ARo sequence probes the effect of nebular geometry through changes in the size of 
the empty central cavity. For AO, AN, and A(N+0), abundances were varied to determine 
the extent to which nebular abundances affect the icp abundance discrepancies. We 
computed ten models, A (log T,,,, log L,,,), in which a "star" moves along a typical Harman- 
Seaton sequence in its PN evolution. That is, a range of stellar UV flux distributions 
(assumed to be blackbodies) and luminosities was considered. For comparison, we also 
considered seven nonLTE models, A(logT*, log which were chosen to approximate 
models from the A(logT^,, logL^) sequence to examine the effects of changes in the shape of 
the input spectrum. These nonLTE models are interpolations of a set of models from Klaus 
Werner (Werner et al. 1991) included in the Cloudy package. 

3. Results 

The model-computed ionization structure and emission line profiles of the baseline 
model are shown in the top and center panels of Figure |^, respectively. The figure is 
characteristic of nearly all of the model nebulae. The resulting abundance discrepancies, 
v4(X, r) , are shown in the lower panel of Fig. |l|. Shaded regions show the areas integrated 
for the low-ionization and neutral regions. Note that the differences between the ionization 
volumes for O"*"^ and Ne"*"^ and and N"*" are cross-hatched in the left and right panels, 
respectively. 

The resultant radial profiles of n'(S"'', r), T'{H^ , r), T'{0^^ , r), and the ionic abundances 
A'{Xiomr) are found to be in such close agreement with their "real" counterparts computed 
by Cloudy that we need not distinguish between them hereafter. This agreement is in 
spite of different atomic constants in the codes used in both methods which only serves to 
reassure us that the errors in the icF method are real and measurable. One exception is 
the Ar"'"^ abundance determined from the 7136A line for which Ionic's results consistently 
differ from Cloudy's computed values by a factor of two. Other exceptions occur, but only 
in extreme circumstances. 



3.1. Global icF abundances measurements. 

The emission line profiles were integrated over each of the four types of regions (entire 
nebula, high-ionization, low-ionization, and neutral) using both volume and line-of-sight 
weighting. These results mimic observations using apertures that isolate and average over 
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the corresponding volumes or lines of sight. From these integrated observables n'(S+), 
T'(N~^), T'(0"^+), A'{Xion), the icF correction factors, the derived abundances A'{X), and 
icF abundance discrepancies A{X) were derived for each zone separately. 

A grid of the line flux ratios integrated over the full nebular volume and corresponding 
to each of the present model sequences is presented in Table |I|. The nebular H/5 luminosity 
and radius for each model (relative to the baseline) are tabulated, as are the densities 
and temperatures that would be derived directly from full-volume integrated emission line 
ratios. 

It is interesting to note how the brightnesses of all of the low-ionization lines 
increase (relative to H/?) along the Harman-Seaton sequence through all ten of the 
A(logT*, logL*) models. By the end of this sequence, the [Nll]/Ha line ratio exceeds 1.5 
(assuming a standard Balmer decrement value of 2.8 for the Hq;/H/3 line ratio). Similarly, 
[Ol]/Ha approaches 25%. These are global integrated values, not just local ones near the 
nebular edge. Comparing these resulting flux ratios to shock models, we note that pure 
photoionization models can give rise to low-ionization line strengths comparable to many 
planar shock models with shock velocities < 200 km s~^ (e.g. Hartigan et al. 1994). Within 
the potential uncertainties of abundances determined with icF methods, the line ratios in 
the low-ionization portion of a planetary nebula can not reliably be distinguished from 
shock models without additional information. 

The H/3 luminosities are relevant for studies of the PN luminosity function. In the case 
that the nebulae are ionization bounded, these luminosities are insensitive to the assumed 
nebular geometry or abundances (see Table |I|). However, as expected, they are excellent 
discriminators of the ultraviolet photon flux of the central star, i.e. its evolutionary state. 

The nebular densities and temperatures in Table |l| exhibit some curious features. 
The [Sll] lines are weighted in favor of regions that are typically only partially ionized. 
Consequently n(S''") underestimates the true nebular electron density (3300 cm~'^ in this 
case) by more than 20%. In addition, T(N"'") can differ from T(0^^) by more than 1000 
°K (~ 10%) due to changes in the nebular temperature profile through the N"*" and O"*"^ 
regions. Accordingly, observed T(N"'") and T{0^^) can be expected to differ in real nebulae 
in which no small-scale temperature fluctuations occur (cf Fig. |I|). The sign of the T(N"'") - 
T(0^+) temperature difference reverses at stellar temperatures around 10^ °K. 

On the one hand, abundance errors are exacerbated unless T(N"'") and T{0^^) are 
used for the computation of low- and high-ionization species, respectively. On the other 
hand, the errors which derive from problems in measuring n'(S"''), T'(N+) and T'(0++) are 
often small compared to other potential problems which are inherent to the icp method. 
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A key assumption in the icF method is that the ionization structure of N and O (and 
hence N"*" and O"*") is identical. However, as seen in Fig. |1|, these ionization zones are 
obviously not cospatial! Although charge transfer locks 0+ to H+, charge transfer is far 
less important in the ionization structure of N"*" (Osterbrock 1989). Hence the ionization 
structures of N+ and 0+ are only weakly coupled. Similarly, the ionization fractions He^, 
O"'""'", and Ne"*"^ are assumed to be roughly equal in the icf method. In the Cloudy models 
these deviate from one another, especially near the edge of the high-ionization region. Here 
too, charge transfer affects O^^ to a greater degree than Ne"*"^. These potential errors in 
the icF method have long been recognized (e.g. Pequignot et al. 1978). 

The substantive issue in this paper is a quantitative analysis of their effect on icF 
abundances, i.e. the magnitude of the discrepancy, A{X). The results of A{X) are 
summarized in Table |^ for the Ai?o and A(logT*, logL*) sequences of models. Results for 
the AO, AN, A(N + O) sequences show only small departures from those of the baseline 
model and are not shown. Similarly, the nonLTE model atmospheres don't change the 
resulting icf abundances to any great degree. In general, the new discrepancies differ from 
the blackbody case by < 5% for unresolved nebulae and < 7% for resolved observations. In 
cases where the blackbody abundances were good (discrepancies < 5%) the nonLTE result 
tends to differ by < 2%. 

Several interesting conclusions can be drawn from inspection of the table. The 
abundances derived for He, N, and O from entire-volume observations (i.e., large-aperture 
observations of ionization bounded nebulae) are generally good. A similar result was found 
by Henry et al. (1996). S and Ar are generally not fitted well (errors > 20%) with Ar 
systematically overestimated by > 20% in volume integrals. Errors for Ne can exceed 50% 
or more and vary greatly with the properties of the central star. 

But the apparent accuracy of the icf method can be deceptive. Note in Fig. |I| the 
degree to which A{X, r) varies with radius. In a volume-average over the entire nebula, 
these variations tend to cancel; however, many of the local values of A{X, r) deviate from 
unity. In regions of low ionizations the icF discrepancies of He, N, Ne, and Ar often 
approach 20%, and sometimes even exceed 100%. Local icF abundance discrepancies are 
considered in 553.2 and 3.3. 



3.2. Radial Icf abundance measurements in the baseline model 

We next consider the abundances derived from observations in which the baseline 
nebula is spatially resolved. Refer to Fig. ^and the baseline model results shown in Table ^ 
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In the outer parts of the nebula where the icF approximation N/0 = N'^/O^ breaks down, 
^(N) increases from 0.7 in the nebular interior to 2 in the neutral zone. If the emission 
lines of the baseline model were observed using a long slit, the use of the icF approximation 
would lead to the erroneous conclusion that both N/0 and N/H increase by a factor of three 
or more! 

The abundances determined solely or largely from lines of high ionization. He and Ne, 
also show large discrepancies in low-ionization regions; however, unlike [Nil], only a small 
fraction of the helium and neon line fluxes arise outside the high-ionization zone. Therefore 
the magnitude of the abundance discrepancy will not be noticeable in practice. Nonetheless, 
icF abundance errors in He and Ne will certainly arise in isolated low-ionization ansae, 
lobes, or halos. 

Similarly, the derived abundances of S and Ar using only the standard [Sll] and [Arm] 
lines are affected. However, their icF discrepancies can be ameliorated if lines of [Slll] and 
[Ariv] are detected and included in the icF abundance analysis. 

Clearly, except for O, all icF abundances are vulnerable to large systematic errors with 
nebular radius, especially in or near the ionization front. This is particularly apparent in 
the lower panel of Fig. |I|. 

3.3. Radial Icf abundance measurements for other models 

The icF discrepancies for other models show some interesting trends from region to 
region. Here we point out only a few of them. Consider first the ARo sequence. As the 
radius of the central hole grows and the impinging stellar radiation becomes more dilute 
(i.e. the ionization parameter decreases), the relative size of the low-ionization zone grows 
and the boundary between the high- and low-ionization regions becomes less distinct. 
This has a profound effect on abundances measured in the high-ionization zone, as seen 
in Table |^. Hollow PNe such as Abell 39 and halos such as NGC 6826, ionized by highly 
dilute stellar radiation, are likely to be the most affected. 

Consider the A(logT*, log L*) sequence. Abundance discrepancies are very sensitive to 
the assumed temperature of the excitation source, as is easily seen in Table |^. The first 
model in the sequence represents a cool stellar ionizing source, logT^, = 4.55. In this case 
the high-ionization zone occupies less than half of the nebular volume, and the boundary 
between it and the low-ionization zone is not sharply defined. In this transition region 
the various ^(X)'s are large. Only oxygen, for which both important ionization stages 
are visible, has a reliably measurable abundance. This serves as a warning for abundance 
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measurements of Hll regions as well as low-ionization PNe such as IC418 and BD+30°3639. 

As the stellar temperature increases in the sequence the ionization patterns seen in 
Fig. are quickly established. The low-ionization zone becomes smaller in radius, and 
the lines that are emitted from it drop in overall importance. As a result, the abundance 
discrepancies A{X) are closest to unity for the second through fifth models in the sequence. 

As the stellar luminosity drops (seventh through tenth models), the low-ionization 
zone again fills more of the nebular volume. Correspondingly, the ^(X)'s increase. As 
discussed earlier, the emission fine ratios of low-ionization species such as 0°, S"*", O"^, N"*", 
and Ar"*"^ resemble values found in low-velocity shocks (Hartigan et al. 1994). At the same 
time, high-ionization lines from He"*", He"'"^, and Ne^^ remain bright. Except for 

oxygen, the ionization correction factors become large and the corresponding abundance 
discrepancies mount. 

4. Discussion 

Ionization correction techniques have long been used in both galactic and extra- 
galactic astronomy. As early as 1985, Bob Rubin recognized problems with derived nitrogen 
abundances. Our work bears out his assertions. We confirm that as electron temperature 
drops the ratio N"'""'"/0"'"^ rises. 

While the present models focus on PNe, the icf abundance discrepancies are important 
for observations of galactic and extragalactic Hll regions (Gonzalez-Delgado et al. 1994, 
Martin 1996, and Storchi-Bergmann et al. 1996 to name only a few recent examples). 
Storchi-Bergmann et al. (1996), in their examination of the central regions of AGNs, 
assume that icF assumptions can lead to an under-estimate in the nitrogen abundance 
of 20%; however, our results suggest that the discrepancy could be 40-50% based on the 
nebular temperatures and densities that they find. 

Nitrogen in radiation-bounded objects like planetary nebulae can be easily 
overestimated, particularly in lines of sight that largely or exclusively sample the low- 
ionization portion of the nebula where the N"*" profile extends beyond both those of H"*" and 
O^. A slit which lies across an ionization-bounded nebula should clearly show a rise in 
N relative to both O and H due to the form of the icp for nitrogen. This effect certainly 
contributes to the results of Balick et al. (1994) and Guerrero et al. (1995), although it 
is unclear whether or not it can completely account for the reported variation in N/0. 
Similar, although less pronounced, problems with Ne and S also exist. Surprisingly, the 
recent results of Corradi et al. (1997) show the expected radial increase in all species except 
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nitrogen. These radial abundance discrepancies will be exaggerated for highly extended 
nebulae such as bipolars in which the flux of radiation striking most of the lobe walls is 
both dilute and incident at high obliquity (see the values for the ARo series in Table 0). 

Since 1967, Peimbert, Torres-Peimbert and their collaborators have been examining 
the physical effects of both large- and small-scale temperature fluctuations on abundances 
determined by the icF method. Recently, Peimbert et al. (1995b) used emission line 
temperatures from different elements to attempt to measure the temperature fluctuations 
in 17 Type I PNe. To quantify the fluctuations, they measured temperatures for numerous 
emission lines and compared them to one another. Our work clearly shows that this is 
not a sound test as the lines originate at different points in the nebula due to different 
sensitivity to temperatures and densities. An examination of the temperatures in Table 3 
of Peimbert et al. (1995b) shows that roughly half of the objects (notably including all 
of the Minkowski objects in their sample) have T(N+) — T(0^+) values that are within 
the differences that one should expect to measure in a smooth, homogeneous PN with no 
small-scale temperature fluctuations. 

Recently, a puzzling factor-of-three-to-flve disagreement in the oxygen abundances 
measured from forbidden and recombination lines has arisen (Liu & Danziger 1993, Peimbert 
et al. 1993, Liu et al. 1995, Mathis 1995). The present study shows that O abundances 
measured using the icf method are generally accurate to ~ 5% measured globally, and 10% 
measured with a pencil beam near the nebular center (before observational errors). Thus 
our results do not resolve this controversy. 

5. Conclusions 

The present model results show that use of icF techniques for abundance analysis for 
spatially unresolved objects gives results accurate to ~ 20% as found by Henry et al. (1996). 
Larger discrepancies can occur for some nebular conditions including cool ionization sources 
or highly dilute incident stellar UV radiation. 

For spatially resolved nebulae, icf abundance errors can be high locally, especially for 
N, Ne, and Ar in or near an ionization front. These errors have led to probable overestimates 
of N/0 abundance variations in FLIERs and bipolar PNe. 

Further, the emission line ratios predicted for the low-ionization species in or near 
ionization fronts in resolved observations resemble those predicted for shocks. This is the 
cause of some ambiguity, if not confusion, in the interpretation of line ratios. 
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The magnitudes of errors in the icF niethod differ for nebulae of various sizes and 
cxcitation/ionization sources. In addition, the icF errors depend on the extent to which 
ionizing photons are absorbed or escape (i.e. whether the nebula is ionization or density 
bounded). Consequently, statistical studies of abundances in many objects need to account 
for these types of conditions. 

Nebulae excited by stars with < 45000°K have relatively large icF discrepancies. 
Our work has concentrated on conditions typical in PNe and thus suggests that a detailed 
study of icF abundances in galactic and extragalactic Hll regions needs to be undertaken in 
order to assess the veracity of icF^derived abundances and gradients in the disks of spiral 
galaxies. 

It is a pleasure to thank the developer of CLOUDY, Gary Ferland. This research was 
sponsored by NSF grant AST 9417112 and NASA GO-6117. 
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A. Empirical Procedures for Deriving Abundances 

A.l. Densities and Temperatures 

Straightforward analytical expressions can be used to derive densities, temperatures 
and ionic abundances from the measured emission line ratios. The functions shown below 
were obtained by a linear least-squares fit to a series of results from the lONiC and Temden 
tasks in the IRAF/STSDAS software package. lA or I(AA) is the measured intensity of a 
line relative to H/? = 1 (not 100) after correction for reddening. 

The task temden in the nebular package in IRAF returns a value of density or 
temperature for a specific set of diagnostic line ratios. The density formula is computed 

assuming T = 10^ °K and is valid for 100 < n(cm^^) < 30000. Assuming that the atomic 
data are correct, the errors of the fits for n(S^) and n(Cl^^) are typically 5% or less. Those 
for n(Ar"'""'""'") are about 8%. As in the paper above, X+ implies X+ZH"*" and X implies X/H. 

76717 
16731 
/5517 
15538 
74711 
74740 



log[n(S+)] = 21.683 - 41.177r + 39.848r^ - 13.954r^ where r 

log[n(Cl++)] = 14.775 - 19.036r + 20.275r^ - 8.3713r3 where r 

log[n(Ar+++)] = 12.713 - 5.7379r - 5.9903r^ + 14.375r^ - 7.1339r^ where r 

The temperature formulae are computed assuming n — lO^cm"^ and are valid for 
6000 < T(°K) < 17000. Errors are less than 1%. 

7(6548 + 6583) 



log[T(N+)] = 5.96r-°-^®^ + 8.13 where r 

log[T(0++)] = 5.22r-°-2°^ + 8.185 where r 



75755 
7(4959 + 5007) 



74363 

Normally some quality-weighted average value of n and T are adopted before proceeding 
further. However, if the data are good, then values of n and T can be used which are 
optimized for the ionization potentials of the ions in the expressions below. 



A. 2. Ionic Abundances 

Ionic abundances of He are computed from expressions fitted to effective recombination 
coefficients found in Osterbrock (1989) (Tables 4.4, 4.5, and 4.6). A density of n = 10^ 
cm~^ has been assumed. We acknowledge that recent work (e.g. Smits 1996) has provided 
better coefficients for helium, however since our goal is to identify errors in the method. 
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rather than differences in atomic coefficients, we have chosen to use Osterbrook's values 
which match those used in CLOUDY. The deviation from curves fit to Osterbrock's tables 
is smaller than 5%. 

The task ionic in the nebular package in IRAF was used to compute various ionic 
abundances from forbidden line ratios (relative to H/5) at T = 7500, 10000, 12500, and 15000 
°K for each of the following densities: 30, 100, 300, 1000, 3000, 5000, 10000, 15000, 20000, 
25000, and 30000 cm~^. Empirical expressions for N, O, Ne, S, and Ar isotopes whose 
coefficients (including those inside the exponentials) are determined from a least-squares fit 
to ionic's results as well as those for He appear in Table ^ These fits are not accurate 
outside of the ranges of temperature and density above. 

In Table |], a; = n/\/T and line ratios are all relative to H/9 = 1, not 100. To the 
right side of the fits, residuals of the fit are shown at 7500, 10000, 12500, and 15000 °K, 
respectively, as a percentage of the fitted value. 

A. 3. Total abundances 

The ionic abundances are multiplied by ionization correction factors to derive total 
abundances. The primary source for the expressions in Table ^ is Kingsburgh & Barlow 
(1994; KB). Alternate sources of ionization correction factors include Koppen et al. (1991; 
KAS), Peimbert et al. (1995b; PLTP) and Freitas-Pacheco et al. (1995; F-PBCI). 
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Fig. 1. — The ionization structure (upper panels; from Cloudy 90.02), line emissivities 
(middle panels) and i^-derived abundance discrepancy (lower panel) for the baseline model. 
The shaded regions are zones where 0°, O"*", 0++ and higher ionization species dominate, 
respectively. The inner edge of the nebula is at 10^^ cm. 

Cross hatching in the upper left panel shows deviations from the icF assumption 
that O"'""'" and Ne^+ are cospatial. Similarly, cross hatching in the upper right panel indicates 
deviations from the assumption that N+ and 0+ are cospatial. In the lower panel, notice the 
tendency for the discrepancies to cancel when integrated over large radii or volumes. This 
implies that icF abundances are more likely to be specious when sampled through smaller 
apertures. 
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Table 1. Line intensities (relative to Hp = 100), Hp fluxes (baseline = 1), and nebular radii (baseline = 1) over entire nebular volume for the various models. 
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3.03 
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2644 
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* The input parameters for baseUne model (shown in boldface) are: logioabundances H:He:C:N:O:Ne:S:Ar=0:-l :-3.52:^:-3.22:-3.82:-4.82:-5; density=0(3000) cm-3 within (beyond) 
Ro=10"cm; central star: logLj=37.4(ergs-i), blackbody spectrum at logT„=5(°K). The resulting baseline values of log L(HP)ergs-i=35.36 and log nebular radius (cm)=17.55. 



Table 2. Abundance discrepancies for the grid of models. 
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1 AA 

I.UU 


1.13 


3.75 


1.00 


1.84 


1.02 


1.00 


1.13 


D. /J 


1 AA 
I.UU 


1 QA 

i.o4 








LOW-IONIZATIONZONE: (x(0+) > 50%) 










1 _ _1 A 'J 

10gKo=lo.3 


1 1 Q 

i.io 


1 OA 
i.ZU 


1.01 


3.95 


1.10 


1.94 


1.17 


1.19 


1 A 1 
1.01 


"2 OA 

j.yu 


1 1 A 
i.iU 


1 CiA 

i.y4 


io. / 


1.17 


i.io 


1.01 


3.59 


1.09 


1.93 


1.16 


1.18 


1.01 


"2 

J.JJ 


1 AO 

i.uy 


1 CiA 

i.y4 


l/.U 


1 1 Q 

l.lcS 


l.ly 


1.01 


3.90 


1.10 


1.94 


1.17 


1.19 


1.01 




1 in 
I.IU 


1 OA 


1 T ■a 
i /.J 


i.lD 


1 n 

1.1 / 


1.01 


3.57 


1.09 


1.94 


1,16 


1.17 


1.01 


1 <o 


1 AO 

I.uy 


1 OA 

i.y4 


17.7 


1.13 


1.11 


1.01 


3.50 


1.09 


1.98 


1,13 


1.11 


1.01 


3.47 


1.09 


1,98 


log(T^,L^) = 


























4.55,37.4 


0.57 


0.80 


0.96 


0.40 


1.34 


1.12 


0.66 


0.75 


0.97 


0.40 


1.28 


1,30 


4.70,37.4 


1.00 


0.82 


0.98 


1.86 


1.12 


2.00 


1.00 


0.81 


0.98 


1.83 


1.11 


2.00 


4.85,37.4 


1.11 


1.06 


1.01 


3.38 


1.08 


1.99 


1.11 


1.06 


1.01 


3.34 


1.08 


1.99 


5.00,37.4 


1.18 


1.19 


1.01 


3.90 


1.10 


1.94 


1.17 


1.19 


1.01 


3.85 


1.10 


1.94 


5.15,37.4 


1.19 


1.22 


1.01 


3.73 


1.10 


1.90 


1.18 


1.22 


1.01 


3.67 


1.10 


1.91 


5.25,37.1 


1.17 


1.21 


1.02 


3.40 


1.09 


1.90 


1,17 


1.21 


1,02 


3.34 


1.08 


1,90 


J.Z.J, JU.U 


1 17 


1 20 


1.03 


3.50 


1.08 


1.91 


1,16 


1.20 


1,03 


3 42 


1 08 


1 91 


< on TA 1 


i.iS 


101 

i.Zi 


1.02 


4.71 


1.11 


1.94 


1.17 


1.20 


1.03 


4. J / 


1 1 A 
i.iU 


1 O/t 

1.94 


D.Uj,jj.u 


1.11 


1 1 

i.iz 


1.01 


4.10 


1.12 


2.03 


1.11 


1.11 


1.01 


1 O"? 

i.y 1 


1 1 

i.iz 


Al 




1 AA 

i.Uu 


1 1 A 
I.IU 


1.04 


5.64 


1.11 


1.97 


1.06 


1.10 


1.04 


j.jU 


1 1 A 
I.IU 


i.yo 










NEUTRAL ZONE*: (x(Oo) > 50%) 










lOgKQ=iO.J 


1 QT 

i.y / 


O OQ 


1.02 


>100 


2.35 


2.64 


1,97 


2.28 


1.02 


1 AA 
>iUU 


O ■25 


z,o4 


io. / 


1 Q'2 

i.y J 


O 01 
Z.Zi 


1.01 


>100 


2.25 


2.60 


1,93 


2.21 


1.01 


^1 AA 
>iUU 


05 

Z.ZJ 


fin 


l/.U 


1 OT 

I.y / 


L.LI 


1.02 


>100 


2.34 


2.64 


1.97 


2.27 


1.02 


>100 





Z.U4 


1 T 1 
i /. J 


1 OO 


1 Q 
Z.io 


1.01 


>1()0 


2.24 


2.58 


1.92 


2.18 


1.01 


^ 1 AA 
> i UU 


0/1 
Z.Z4 


Z.JO 


17.7 


1.90 


2.03 


1.01 


>100 


2.25 


2.61 


1.90 


2.03 


1.01 


>100 


2.25 


2.61 


log{T*,L*) = 


























4.55,37.4 


0.00 


2.27 


1.14 


>100 


2.27 


0.00 


0.00 


2.27 


1.14 


>100 


2.27 


0.00 


4.70,37.4 


1.77 


2.48 


1.06 


>100 


2.39 


2.67 


1.77 


2.48 


1.06 


>100 


2.39 


2.67 


4.85,37.4 


1.83 


2.19 


1.02 


>100 


2.11 


2.58 


1.83 


2.19 


1.02 


>100 


2.11 


2.58 


5.00,37.4 


1.97 


2.27 


1.02 


>100 


2.34 


2.64 


1,97 


2.27 


1,02 


>100 


2.34 


2.64 


5.15,37.4 


1.93 


2.12 


1.01 


>100 


2.25 


2.60 


1.93 


2.12 


1.01 


>100 


2.25 


2.60 


5.25,37.1 


1.92 


2.05 


1.01 


>100 


2.25 


2.59 


1.92 


2.05 


1.01 


>100 


2.25 


2.59 


5.25,36.6 


1.95 


2.02 


1.01 


>ro() 


2.31 


2.65 


1.95 


2.02 


1.01 


>100 


2.31 


2.65 


5.20,36.1 


1.91 


1.93 


1.01 


>100 


2.24 


2.60 


1.91 


1.93 


1.01 


>10() 


2.24 


2.60 


5.05,35.6 


1.92 


1.91 


1.01 


>100 


2.33 


2.69 


1.92 


1.91 


1.01 


>100 


2.33 


2.69 


5.00,35.1 


1.85 


1.80 


1.01 


>100 


2.26 


2.65 


1.85 


1.80 


1.01 


>100 


2.26 


2.65 



























* Underlined entries indicate that the corresponding lines are very weak or confused by much brighter 

nearby regions, 
t Errors are approximately differences of adjacent rows, 
t Model truncated at x(H+) = 30%. The neutral zone is notcomplete. 



Table 3. Ionic abundances relative to 





Fitted Expression 


Fractional error at 
various temperatures* 


Notes 


He*/H* = 


I(5876)/I(H|5) 0.739 ( T/10,000 f™ 


<2% 






He*/H* = 


I(6678)/I(HP) 2.578 ( T/10,000 f^" 


<2% 






He'"/H* = 


: I(4686)/I(HP) 0.0832 ( T/10,000 )° ^'' 


<2% 






n*/h'' = 


I(6583)/I(HP) 1.341 x 10"' exp (26,485/1) / (1 + 1.2224x 10"' x) 


2% 2% 1% 


2% 




0°/H* = 


I(6300)/I(H|i) 4.80x 10"* exp (37,076/T) * (0.383 + 1.2 x 10"* T - 5.443 x 10"' T^) 


1% 1% 1% 


1% 


1 


0* / H* = 


I(3726+3729)/I(HP) 4.684« 10"' exp (43,706/T) * (1. 14 - 1.156x 10"' T) * (1.05 + 1.534x 10"" n - 2.54x 10"'V) 


1% 1% 1% 


1% 




0*/H* = 


I(7325)/I(HP) 4.20x 10"' exp (62,545/1) * (0.612 + 0.405 [e ' - ( 1.1 - 0.9x 10"'t) e"]/[e"' + e*']) , z = (logion - 3.04) 


7% 2% 2% 


2% 


2 


0**/H* = 


= I(5007)/I(H|J) 1.152x10"' exp (34,609/T) /(I + 1.2487 xlO"^) 


4% 2% 2% 


3% 






= I(3869)/I(HP) 1.33x10"* exp (41,250/1) 


3% 2% 1% 


1% 




s*/h' = 


I(6717+6731)/I{HP) 1.660x 10"' exp (25,079/1) * (1.03 + 1.89x lo"*!!) 


3% 1% 1% 


2% 






1(63 12)/I(H3) 3.1x10"' exp (43,250/1) / (1 + 6.93 x 10" V") 


2% 2% 2% 


2% 




At* I H* = 


= I(7136)/I(HP) 8.189X 10"' exp (24,243/T) 


2% 2% 1% 


2% 






= I(4740)/I(HP) 7.433 x 10"' exp (29,789/T) * (1 + 0.47e"" ""'*) 


1% 1% 1% 


1% 




4+ + 

Ar /H = 


= I(7006)/I(H|3) 1.325x10"' exp (26,279/1) *{ 1.005+ [1.15x10"^ n - 0.049] exp f -0.115 n""'' )) 


6% 1% 1% 


1% 





* The values of density n and temperature T used in these expressions for ionic abundances should be derived from lines of ions with similar ionization 
characteristics, x = n/ifj- Errors in the second column are computed at 7500, 10000, 12500, and 15000 °K, respectively, as described in the text. 
Notes to Table 3: 

1. The fit is good; however, IONIC is inaccurate owing to the complex form of the target area parameter CI. 

2. The 7325A line of consists of a blend of four lines at 7321.7, 7322.3, 7332.2 & 7332.8 A. 



Table 4. Total abundances relative to H when only optical Unes are measured* 



icF Expression Ref. 

He = He+ + He++ 

O = (O* + O**) (He / He*)™ KB 

N = ]Sr(0/0*) KB 

Ne = Ne"(0/0") KB 

S = (S'' + S**)[l-(l-0*/0)']"' KB 
where, if S++ is not measurable, adopt S++ = S+ [4.677 + (0++/0+)» «3 ] 

S = (S* + S**) [1.43 - 0.196 (0**/0*)'-^' ] KAS 

Ar=1.87Ar** if the only measurable line is [Ar Till KB 

Ar = 1 .34 Ar** ( 1 + 0*/0**) if the only measurable line is [Ar III] F-PBCI 

Ar = (Ar** + Ar'* + Ar*^(l-N+/N)-i PLT-P 

* In the present paper only the expressions of KB have been adopted. 
The [S in]6312A is assumed to be measurable. 



